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*T3 ' Abstract 

S3 ■ 

An isotropic elastic half space is prestrained so that two of the 

principal axes of strain lie in the bounding plane, which itself remains 

free of traction. The material is subject to an isotropic constraint 

of arbitrary nature. A surface wave is propagated sinusoidally along 

the bounding surface in the direction of a principal axis of strain and 

JO | decays away from the surface. The exact secular equation is derived 

by a direct method for such a principal surface wave; it is cubic in a 

quantity whose square is linearly related to the squared wave speed. 
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i^) For the prestrained material, replacing the squared wave speed by 

C*~) ' zero gives an explicit bifurcation, or stability, criterion. Conditions on 

the existence and uniqueness of surface waves are given. The bifurca- 
tion criterion is derived for specific strain energies in the case of four 
S^ | isotropic constraints: those of incompressibility, Bell, constant area, 

H ■ and Ericksen. In each case investigated, the bifurcation criterion is 

found to be of a universal nature in that it depends only on the prin- 
cipal stretches, not on the material constants. Some results related to 
the surface stability of arterial wall mechanics are also presented. 



1 Introduction 

There has been a large body of literature on surface waves propagating si- 
nusoidally along the bounding planar surface of an elastic half space with 
attenuation of the wave amplitude in the direction normal to the bounding 
plane. The earliest paper is that of Rayleigh [T], who investigated surface 
waves propagating across the surface of an isotropic earth in the context 
of seismology. The modern theory of surface waves derives in large mea- 
sure from the sextic formalism of Stroh [2]. This approach has been em- 
ployed by many different authors to address many different problems of linear 
anisotropic elasticity and the results have been comprehensively reviewed by 
Ting [3]. In parallel, tremendous progress has been made in the theory of 
small-amplitude surface waves propagating on finitely deformed, nonlinearly 
elastic half-spaces, from the seminal works of Hayes and Rivlin [3] (compress- 
ible materials) to those of Dowaikh and Ogden [5] and many others. 

In the present paper we consider an isotropic elastic half space prestrained 
so that two of the principal axes of strain lie in the bounding plane which 
itself remains free of traction. Additionally, the material is subject to an 
isotropic constraint of an arbitrary nature. One example of such a constraint 
is incompressibility, which has received much attention in the literature, and 
another is the Bell constraint, whose experimental and theoretical properties 
have been thoroughly reviewed by Beatty [6, Chap. 2]. A wave is propagated 
sinusoidally in the direction of one of the principal axes lying in the bounding 
plane, has no in-plane displacement component in the direction of the second 
principal axis lying in this plane, but has attenuating amplitude in the di- 
rection of the third principal axis, orthogonal to this plane (principal surface 
wave). The theoretical framework for the study of these waves is set up in 
Section 2, including a discussion of four example of isotropic constraints, the 
incremental equations of motion and of constraint, and the strong ellipticity 
condition. 

In Section 3 we derive an explicit secular equation for these principal 
surface waves that is cubic in a quantity whose square is linearly related to 
the squared wave speed. On restricting attention to an unstrained isotropic 
material we find that the secular equation reduces to that found by Rayleigh 
[1] in the incompressible case. We show that for the unstrained material all 
isotropic constraints are the same, that is all reduce to incompressibility [7]. 
Returning to the general secular equation of the prestrained material and 
replacing the squared wave speed by zero we obtain an explicit bifurcation, 



or stability, criterion for the material. Some results on the existence and 
uniqueness of surface waves are given. 

Finally, in Section 4 we examine the bifurcation criterion for each of 
four examples of isotropic constraint and obtain explicit results by choosing 
special forms of the strain energy function. In all cases considered we find 
that the bifurcation criterion is of a universal nature in that it depends only 
on the principal stretches, not on the material constants. 

2 Preliminaries 

2.1 Deformed constrained half-space with a free plane 
surface 

We consider a semi-infinite body made of homogeneous isotropic hyperelastic 
material at rest in a configuration B u with strain energy density W per unit 
volume of B u and mass density p. Let (O, X 1? X 2 , X 3 ) be a fixed rectangular 
Cartesian coordinate system such that the body occupies the region X 2 > 0. 
The orthonormal set of vectors {i, j, k} is aligned with the coordinate axes. 
Loads Pi, P 2 , P3 are applied at infinity to deform and maintain the half- 
space in a static state B e of finite pure homogeneous deformation, with corre- 
sponding stretch ratios Ai, A2, A3 in the i,j,k directions. Thus, the position 
of a particle at (X±, X 2 , X 3 ) in B u is at (xi,x 2 ,xs) in B e , where x\ = A1X1, 
x 2 = \ 2 X 2 , x 3 = A3X3. The constant deformation gradient associated with 
the deformation is 

F = Axi^i + Aaj^j + Aak^k. (2.1) 

In an isotropic hyperelastic material the strain energy is a symmetric func- 
tion W(Xi, X 2 , A 3 ) of the principal stretches, i.e. its value is left unchanged 
by any permutation of the stretches Ai, A2, A3. The material is subject to an 
isotropic internal constraint, written as 

r(A 1 ,A 2 ,A 3 ) = 0, (2.2) 

in which T is a symmetric function of the principal stretches Aj. We restrict 
attention to constraints such that: 

Ti > 0, where T; := dT/dXi. (2.3) 



Four examples of such constraints are treated explicitly in this paper and 
they are henceforward denoted by roman numerals: the incompressibility (I), 
Bell (II), areal (III), and Ericksen (IV) constraints, 



F 1 := AxA 2 A 3 -1 = 0, 






T 11 : 


= Ax + A 2 + A 3 - 3 = 0, 


r := A1A2 + A2A3 + A3A1 - 


-3 = 


= 0, 


p/V 


:= A? + Al + A§-3 = 0. 



(2.4) 

The incompressibility constraint is often used for the modelling of finite de- 
formations of rubber-like materials and shows good correlation with experi- 
ment (see for instance Ogden [SJ Chap. 7]). The Bell constraint was found to 
hold experimentally over countless trials on polycrystalline annealed solids, 
including aluminum, brass, copper, and mild steel, see Beatty [6] Chap. 2]. 
The areal (or constant area) constraint has the interpretation that a mate- 
rial cube in the reference configuration with edges parallel to the principal 
axes of strain retains the same total surface area after deformation; it was 
studied from a purely mathematical point of view by Bosi and Salvatori [9]. 
Finally, the fourth constraint was proposed by Ericksen [10] to model the 
behaviour of certain twinned elastic crystals. Although Ericksen proposed a 
multi-constrained model, pursued by Scott [11] in the context of wave propa- 
gation, some authors [T2J [T3] refer to the single constraint ( 12. 4p ^ in nonlinear 
elasticity theory as 'Ericksen's constraint'. 

The general constraint ( 12. 2 p generates the workless reaction tensor N, 
see [HJ, given by N = J^ 1 F(<9r/<9F) T , where J = AiA 2 A 3 . Explicitly, the 
non-zero components of N are 

Nu = J~ 1 \ l T i (no sum). (2.5) 

Thus for the four examples of constraints (12. 4p we find the following con- 
straint tensors: 



N 7 = l, N IJ = J- 1 V, N m = ./^[(trVyV-V 2 ], N IV = 2J~ 1 V 2 , 

(2.6) 
in terms of the left stretch tensor V = diag (Ai, A2, A3). 

Associated with the deformation B u — > B e is the Cauchy stress tensor a 
which takes diagonal form with non-zero components 

a ii = J- 1 \ i W i + PN ti (no sum), (2.7) 

where Wi := dW/dXi and P is a scalar to be determined from the equations 
of equilibrium and boundary conditions as follows. First, we note that for 



P constant, the equations of equilibrium dWij/dxj = are automatically 
satisfied. Next, we assume that the surface x 2 = is free of tractions; it 
follows that a 22 = (and P 2 = 0) and so 

P = -J- 1 \ 2 W 2 /N 22 = -W 2 /T 2 . (2.8) 

Consequently, the constant loads Pi, P3 needed at infinity in order to main- 
tain the half-space in the deformed configuration B e are (k — 1, 3 no sum), 

P k = ~a kk = (W 2 /T 2 )N kk - J~ l \ k W k = (W 2 T k - W k T 2 ) X k /(JT 2 ). (2.9) 

Note that when the boundary surface is not free of tractions (a 22 7^ 0) 
then P in ( I2.8P must be replaced by P = (Ja 22 — \ 2 W 2 )/(X 2 T 2 ). However, 
the primary interest of this paper is in studying the influence of internal 
constraints upon the propagation of surface waves rather than the influence 
of pre-loading, and we assume henceforward that (12. 8p holds. 

2.2 Superposed infinitesimal motion and strong ellip- 
ticity condition 

We now consider the propagation of an incremental motion in the deformed 
half-space B e — )■ B t , described by 

x = x + eu(x,t), (2.10) 

where x is the position in Bt of a particle which was at x in B e and u is 
referred to as the displacement vector. The parameter e is small, so that 
terms of order higher than one in e may be neglected. The linear fourth- 
order instantaneous elasticity tensor B* associated with the motion is [TJJ 
(2.19)] 

Bijki = Bijki + PBij k i = Bij k i — (W 2 /T 2 )Bij k i, (2-11) 

where the non-zero components of B and B are given by [SJ (6.3.15)] 

T „ AjW j — AjWa . n T ~ A;l ,■ — Ail 1 , o , . 

JB m = A ;_ A 2 X l JB ™= A 2_ A 2 Xl {t ^ j) (212) 

ijji " ijij AiVVi) J JDijji J -Dijij Ajl j, (I y^ J ) 

and there is no summation over i or j. Here the second line holds when 
* 7^ Ji A i 7^ Aj] i n the case where i ^ j, Aj = A,,-, it must be replaced by 



" -*-*ijij 2 v iiii Hjj ' "i'*i)'i " ijij ~ 2 V i*ii iijj "^ "i*-i) , > See 

P, (6.3.16)]. Because of hyperelasticity the symmetries B*^ kl = B kli ^ hold so 
that in particular we have 

B UJ3 = B JJH> B hi = B jHJ ( UO SUm )- (2.13) 

The incremental nominal stress associated with this motion is s given by 
[11(3.10)] _ _ 

Sij = on + B* jkl ui jk + pN ij} (2.14) 

where p represents the increment in P. 

The equations of motion, together with the incremental constraint, read, 
from [14, (3.13)], 

~pUj,tt = s ijti , iVii«i,i + N 2 2U 2 ,2 + N 33 u 3 , 3 = 0. (2.15) 

For future convenience, we assume that B* is strongly elliptic, that is, 

B*jki m i n j m kni > 0, for all m, n such that rriiNijnj = 0. (2.16) 

The vectors 

m = iVV cos 9 i + N~j sin j = J^ [(AiTi) - ^ cos 9 i + (A 2 r 2 )^ sin 9 j], 

n = -iV^sin0i + iV^cos0j = ^[-(XiTx)'^ sin0i + (A 2 r 2 )^ cos0j], 

(2-17) 
< < 2ir, are two vectors satisfying ( I2.16P 9. Introducing the quantities A, 
B, C, defined by 

a = (AiA 2 r!r 2 )~ B{ 212 , c = (AiA 2 r!r 2 ) _ 5 2121 , (2.18) 
B — 2 [(^i r i) _ ^1111 + (A2r 2 ) _ 5 2222 ] - (AiA 2 rir 2 ) _ (-Bi 122 + b* 1221 ), 

we obtain from (I2.16P the inequality 

A cos 4 6 + 2B sin 2 cos 2 9 + C sin 4 > 0, (2.19) 

holding for all 0. In particular, the choices 9 = 0, ir/2, arctan(— B/C)^ 1 ^ 2 , 
give in turn 

A>0, C>0, B + y/AC>0. (2.20) 

Using different notations, similar inequalities are given by Knowles and Stern- 
berg [15] for unconstrained materials, by Ogden [8] for incompressible mate- 
rials, and by Destrade [16] for Bell materials. 



3 Surface waves 

Here we specialize the equations of motion to the consideration of an inho- 
mogeneous principal plane wave traveling over the free surface x~2 = and 
derive the corresponding secular equation. 

3.1 Equations of motion and boundary conditions 

A surface (Rayleigh) wave travels sinusoidally with time t in a direction par- 
allel to the plane x~2 = (X 2 = 0) leaving it free of tractions and decaying 
away from this surface as £2 — > 00. For simplicity, we consider a wave propa- 
gating with speed v and wave number k in the principal direction of prestrain 
OXi = Ox\. We refer to this as a principal wave. For this wave, antiplane 
strain decouples from inplane strain, and the displacement components are 
of the form: 

u . = Ui ( k x 2 y^i-vt) ( z = i ?2 ), u 3 = 0, p = kQ(kx 2 )e lk ^- vt \ (3.1) 

in which U\, U2 and Q are functions of k~x~2 to be determined. Similarly, the 
antiplane stress decouples from the plane stress, and the equations of motion 
(12. 151) reduce to 

Sll.l + S 2 1,2 = pUl,tt, Sl 2) l + S 2 2,2 = pU2,tt, N U U 1A + N 22«2,2 = 0. (3.2) 

Because the Wij terms in the components (12.141) of Sij are constant, we 
write the relevant stress components in the form 

s ij = <T ij + kS ij (kx 2 )e ik & 1 - v£ > (i,i = 1,2), (3.3) 

and the equations of motion reduce further to 

iS u + S' 21 = -/w 2 £A, iS l2 + S' 22 = -pv 2 U 2 , zAiIW! + A 2 r 2 ^ = 0. (3.4) 

Explicitly, the S^ are given from f)2.14|) and (I3.3P by 

Sii = iBlmU-L + 5ii22^2 + QN\x, S X2 = B* l22l U[ + %B\ 2X2 V 2 , 

S 22 = iB\ l22 U x + B* 2222 U' 2 + QiV22, S 2 i = B* 2121 U[ + iB* 221 U 2 , (3.5) 

in which the symmetries H2.13f) have been used. 



Now we use (j3.4p and (13.51) to formulate the problem as a system of four 
first-order differential equations for the unknown functions Ui, U 2 , S 2 i, S 22 : 



U[ = -i 



B*. 



B. 



1221 U 2 



2121 



R* 
- D 2121 



S- 



21) 



^ 



.Axr 



s' 



21 



■°1111 



Air 



A 9 r 



2 1 2 



-£* 



122 



AiTi 

A2r2 



A2r2 



R* 

-°2222 



'-Uu 



pv 



On 



22 



R* R* 

- D 2121- D 1212 



R* 2 
- D 1221 



R* 

-°2121 



pf 



p* 
,-D 10 



U 2 _ j^i 5 . 

R* 



21- 



''2121 



U-i 



.Air, 

^2-L 2 



(3.6) 



Equations (j3.6P i and ( I3.6j) q are obtained from ( 13. 5p ^ and ( I3.4p ? . respectively, 
and ( 12. 5 p is employed in the latter. Equation f!3.6|) ^ is obtained by eliminating 
Q between A3.5J) i and (I3.5p ^ using ( 12. 5p . and then eliminating S n between this 
equation and ( I3.4l) i. Finally ( I3.6p ,i is obtained by eliminating S\ 2 between 
( 13.41) 9 and (l3.5[) o and then using ( I3.5lh to eliminate U[ in favor of S 2 i. 

Equations (13. 6p are subject to the boundary conditions of decay as x 2 — > 
00 and of vanishing traction on x 2 = 0: 



S 21 (0) = S 22 (0)=0. 



(3.7) 



Furthermore, each unknown may be written in terms of a single unknown 
function <p of the variable 

z = ——kx 2) 

^2-L 2 

with prime now denoting differentiation with respect to z: 



U! = i(p'(z), U 2 



S> 



22 



+ 



A2r2 



<p{z), S21 
B 2121 <p'"(z) 



^±B* 212l¥ >"(z) + iB{ 221 ip(z), (3.8) 



Aor 



2 j- 2 



Aifj 



p* 
-Dim 



Air 



A 2 r 2 



-(-ST901 + 25 



1122, 
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AiTi 

A 2 r 2 



B, 



2222 
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^W 



Here, the expressions for C/i, £21, and S 22 were obtained from (J3.6P ? 1 > 3, re- 
spectively. The last equation of this set, namely (!3.6P /i. then yields a differ- 
ential equation for <p(z): 



7V"" - (2/3* - pw V + (a* - pV)<p = 0, 



(3.9) 



Wi\\. r? 



where 

" ~~ - D 1212' / _ I \ p I " D 2121 — p2" ' 

a r a p \ 2 

2/3* = £?* 1U - 2-^-^(^22! + S* 122 ) + f ^r J ^2222- (3-10) 

We have used the property X 2 B* i:Ji = \ 2 B*-^ (i ^ j, no sum) derived from 
f !2.12p . By comparison with the quantities A, B, C, defined in ( I2.18p . and 
use of the consequences (j2.20p of the strong ellipticity (S-E) condition (I2.16p . 
we find that these coefficients satisfy the inequalities 



a*>0, 7*>0, (3* + ^a^f>0. (3.11) 

3.2 Secular equation and bifurcation criterion 

We now derive the exact form of the secular equation for a surface wave 
traveling in a principal direction of a deformed isotropic material subject to 
a single isotropic constraint. Because the wave amplitude decays as z — > oo 
away from the free plane z = 0, we seek a solution for cp in the form 

<p(z) = A 1 e- S ' z + A 2 e~ S2Z , &(*) > 0, Sl ^ s 2 . (3.12) 

From (13.91) , the s , are roots of the biquadratic 

7 V - (2/3* - pv 2 )s 2 + (a* - -pv 2 ) = 0. (3.13) 

The roots s 2 of this real quadratic are either both real (and, if so, both 
positive because we must have 9ft(sj) > 0) or they are a complex conjugate 
pair. In either case, sfs 2 , > and so, by (I3.13p . 

0<-pv 2 <a*. (3.14) 

Although necessary, this inequality on the squared wave speed is not sufficient 

to ensure the decay of the wave amplitude, as is seen in the next subsection. 

The boundary conditions ( 13. 7p yield, using (j3.8j) a a. ( I3.10p . and (I3.13p . 

(l*4 + ^M*22l) M + ^*S 2 2 + ^M* 221 ) A 2 = 0, 

si (l*4 + ^M*22l) A t + , 2 (i>4 + ^^1221) A 2 = 0. (3.15) 



The vanishing of the determinant of this homogeneous system gives the sec- 
ular equation, that is, the equation for the wave speed. Introducing the 
quantities 

a* - pv 2 AiTi A 2 r 2 t 

V =\ ; , o = ——B 1221 = ——-f , 3.16 

V 7 A 2 r 2 Aihi 

and using (I3.13P we obtain from the vanishing of the determinant of (I3.15p . 

28* 4- 28* — n* fr* 2 

f(r,*) := r?* 3 + rf 2 + ^ +2 ; % * - 6 — = 0, (3.17) 

after removing the s% — s 2 factor. Equation ( 13. 17ft is the required secular 
equation. 

In the special case where the biquadratic (13. 13[) has double roots, so that 
s 2 , = s 2 , the above derivation is no longer valid. The form (13. 12[) of solution 
must be replaced by 

<p(z) = (A, + A 2 z)e- S ' z 

but it can still be shown that the secular equation is given by A3 . 1 7|) . Thus 
( I3.17P furnishes the secular equation in all cases. 

The cubic equation ( 13 . 1 7[) in r]* has remarkable features: the coefficients 
of the two highest powers are each always equal to unity, irrespective of the 
pre-deformation, constraint, and strain energy function. By ( 13. 161) ^. the term 
independent of rf depends only on Ai, A 2 , and the constraint, but not on the 
strain energy function. For instance, for the I- IV constraints ( \2A\i , we find 
that 

— =1 — — ^ ^ 3 ^ 2 ^ /O -|0\ 

7* ' Ai' A 2 + A 3 Ai' XV { } 

respectively, whatever W may be. From ( 13.10^ 9 it is clear also that «*/7* 
is independent of W and it follows that the secular equation (I3.17P depends 
on W only through the term /3*/7* appearing in the coefficient of the term 
linear in 77*. 

In an undeformed material (Ai = A 2 = A3 = 1), we have a* = /3* = 7* = 
5* = n, where \x is the infinitesimal shear modulus, and the secular equation 
( I3.17P reduces to 



f( v ) := rf + rf + 3r) - 1 = 0, rj = y/l - pv 2 /p. (3.19) 

The unique positive real root of this equation corresponds to pv 2 /p ~ 0.9126, 
in accordance with Rayleigh's result [1] for incompressible linear isotropic 
elastic materials. 
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Equation (I3.19P is valid for all isotropic constraints provided that the 
isotropic elastic material is undeformed. In fact, in the undeformed state 
we can show that all isotropic constraints are equivalent by arguing as fol- 
lows (see Podio-Guidugli and Vianello [7] for an alternative treatment). The 
constraint is T(Ai, A2, A3) = where the function T is symmetric in its argu- 
ments, i.e. its value is left unchanged by any permutation of the stretches 
Ai, A 2 , A3. For small strains write Aj = 1 + e^ where e, is the extension ra- 
tio. The constraint holds for Aj = 1 and for Aj = 1 + e, (for each i) so an 
application of Taylor's theorem gives 

dT dT dT 

-ei + 7T^e 2 + 7T^e 3 = 0, 



<9Ai d\ 2 <9A S 



where terms quadratic in e, are neglected. The partial derivatives are evalu- 
ated at Aj = 1 and, because of the symmetry condition on T, are all equal. 
Thus, in the undeformed state, each isotropic constraint takes the form 

ei + e 2 + e 3 = (3.20) 

of the constraint of incompressibility for infinitesimal deformations. 

The bifurcation criterion is obtained by writing v = in the secular 
equation (13.171) and indicates when the half-space might become unstable 



(3.21) 



<m- 


/a*\5 a* 2(3* + 25* -a* 


la* 


6* 2 

/y*2 


= 


which becomes 


r| n /3* + 5*T 2 5* 2 

yj + 2 7 * r\ 7* 2 "" 










(3.22) 

on using fl3T0l) o to eliminate a*/j*. Using (pTTjl o. f l2~T2l) . fEUOl) and (13T6]U . 
we may rewrite (13.221) in terms of the derivatives with respect to Aj of the 
strain energy function W and the constraint T: 

T 2 2 W n - 2T 1 T 2 W 12 + T 2 W 22 
+ T l {T 2 W 1 - T 1 W 2 )/X 1 - (T n T 2 2 - 2T 1 T 2 T 12 + T 22 T\)W 2 /T 2 = 0, (3.23) 

an explicit form of the bifurcation criterion. 
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3.3 Existence and uniqueness of a surface wave 

Following Chadwick [18] we recast our secular equation (I3.17P as 



f(v) = det 



5* 
5* 

7]* 



r 



where 



{r}* + iy 



a* + 7* - 2/3* 

ry* 



. (3.24) 



Chadwick used this "matrix reformulation" of the secular equation in the case 
of incompressible (I) materials (where 5*/7* = 1) to derive, in a rigorous 
manner, essential results about the existence and uniqueness of a surface 
wave. For a general isotropic constraint ( 12. 2p . <5*/7* is n °t necessarily equal to 
1, but his method, together with the S-E inequalities ( 13. lip , can nevertheless 
be directly transposed to our problem. In summary, we obtain the following 
results. 

The limiting speed v, which is the upper bound of the subsonic interval 
I = [0, v] for v where the decay requirement 9£(sj) > is satisfied, is defined 
by 



pv 



a 



when 2/3* > a*, 
2[/3* - 7* + v/7*(a*-2/3* + 7*)] < a*, when 2/3* < a*. 



(3.25) 



A necessary and sufficient condition of existence for a root vr in I is 

/(0)>0, f(v)<0. (3.26) 

Finally, when a root exists, it is unique. 



4 The bifurcation criterion in special cases 



We specialize the bifurcation criterion (I3.23P to each of the four examples of 
isotropic constraints considered in this paper and pick special forms of the 
strain energy in order to obtain explicit results. In each case we find that 
the bifurcation criterion is of a (relative) universal nature in that it depends 
only on the principal stretches, not on the material constants. 
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Incompressibility (I) 

For incompressibility, the bifurcation criterion (I3.23J) reduces to [5]: 

\{W XX - 2\ 1 \ 2 W 12 + X 2 2 W 22 + X 2 [W X + (2 - X[ X X 2 )W 2 \ = 0. (4.1) 

We consider the strain energy function for "generalized Varga materials" , 

W = di(Ai + A 2 + A 3 - 3) + d 2 (AxA 2 + A 2 A 3 + A 3 Ax - 3), (4.2) 

in which d\ and d 2 are material constants. The S-E conditions (13.111) lead to 
d\ + 2<i 2 A 3 > 0, and so d x > 0, d 2 > or d\ > 0, d 2 > 0. The bifurcation 
criterion (14.11) leads to 



(d 1 + \ 3 d 2 )(3\ 1 -\ 2 ) = 0, or A 2 = 3Ai, (4.3) 

independently of the choice of material constants. 

Bell's constraint (II) 

For a Bell constrained material, (14. ip reduces to the bifurcation criterion [T5] : 

Ai(Wn - 2W 12 + W 22 ) + W^ - W 2 = 0. (4.4) 

For the specific example of the "simple hyperelastic Bell material" [6] we take 

W = d 2 (AiA 2 + A 2 A 3 + A 3 Ax - 3) + d 3 (A!A 2 A 3 - 1). (4.5) 

Here the S-E conditions ( 13 .lip lead to — (d 2 + 2<i 3 A 3 ) > 0, and so d 2 < 0, 
d 3 < ot d 2 < 0, d 3 < 0, while the bifurcation criterion (14. ip reduces to 

(4 + A 3 4)(A 2 -3A 1 )=0, (4.6) 

leading again to (J4.3P . 

Areal constraint (III) 

The bifurcation criterion (14.11) reduces to 



(Ax + A 3 ) 2 iy n - 2(Ai + A 3 )(A 2 + X 3 )W 12 + (A 2 + A 3 ) 2 iy 22 
+ A^ X (A X + A 3 )(A 2 + A 3 )Wi + (A 2 + A 3 )[2 - X^(X 2 + X 3 )]W 2 = 0. (4.7) 
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For the following material, which we term "simple hyperelastic areal mate- 
rial", 

W = di(Ai + A 2 + A 3 - 3) + 4(AiA 2 A 3 - 1), (4.8) 

the S-E conditions (13.111) lead to d\ — d^X^ > 0, and so d\ > 0, d 3 < or 
d\ > 0, d 3 < 0. The bifurcation criterion (J4.7P reduces to 

(di-A§d 3 )(3Ai-A 2 ) = 0, (4.9) 

leading once more to (14.31) . 

Ericksen's constraint (IV) 

For Ericksen materials the bifurcation criterion (14. ip reduces to 

\\W n - 2\ 1 \ 2 W l2 + X\W 22 + A 2 [W r i - AJ 2 (A? + AiA 2 + \ 2 2 )W 2 ] = 0. (4.10) 

For the following "simple hyperelastic Ericksen material" , 

W = D 2 {\\\1 + A^ + >^\l - 3) + D z {\\\l\l - 1), (4.11) 

the S-E conditions (I3TTT]) lead to -(D 2 + D 3 \l) > 0, and so D 2 < 0, D 3 < 
or D 2 < 0, D 3 > 0. The bifurcation criterion (I4.10p simplifies to 

\l + 5Xl\ 2 - X^l - Xj = 0. (4.12) 

Summary 

For the four specific constraints and strain-energy functions above, we found 
relative-universal bifurcation criteria, even though two independent material 
constants were involved. Moreover, the bifurcation criterion for each of the 
first three forms of W turns out to be the same, although the constraints are 
different; of course, the corresponding critical stretch ratios differ in each case 
because they are obtained by solving the bifurcation criterion in conjunction 
with the constraint condition. These differences are highlighted in Table 1, 
where the critical stretch ratios (Ai) cr for surface stability in compression of 
the four constrained materials presented in this Section have been computed 
in the case of plane strain A3 = 1, and of equibiaxial strains X 2 = A3 and 
Ai = A3. All in all, it seems that simple hyperelastic Bell materials can 
be compressed the most before the bifurcation criterion is reached, whilst 
Ericksen materials of type ( 14. lip can be compressed the least. 
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Table 1: Critical stretch ratios (Ai) cr for surface stability 

A2 = A3 A3 = 1 Ai = A3 



simple Bell 


0.429 


0.500 


0.600 


simple areal 


0.447 


0.535 


0.655 


generalized Varga 


0.481 


0.577 


0.693 


simple Ericksen 


0.603 


0.658 


0.730 



Some questions worth investigating are raised by the previous results: 
what is the largest class of energy function for which the result is univer- 
sal for any given constraint? for what energy functions is the bifurcation 
criterion the same for two (or more) different constraints (separately or to- 
gether)? for a material whose energy function is linear in the invariants of 
the stretch tensor such as (14.21) . (14. 5p . (14. 8p . which constraints lead to the 
same bifurcation criterion ( 14.3J) ? and so on. Unfortunately we must, because 
of space limitations, leave these problems open. 

Concluding remark: an incompressible model for human thoracic 
aorta 

We conclude with an example of a non-universal bifurcation criterion using a 
model taken from the biomechanics literature. Horgan and Saccomandi [20] 
recently discussed the merits of the "limiting chain extensibility model" pro- 
posed by Gent [2T] and its applicability to the modelling of strain-stiffening 
biological tissues. For this incompressible (I) material, the strain energy is 

W = ^J m ]nfl- Xl + X2 ^ Xl ~ 3 y \j + \l + \l<3 + J m , (4.13) 

where /1 is the shear modulus and J m is a constant. The smaller J m is, the 
stiffer the material becomes; conversely the body becomes more deformable 
as J m increases, with the neo-Hookean model as a limit for J m — > 00, where 
there are no restrictions on the values that the stretch ratios may take. 

For the Gent model (I4.13P the S-E conditions (13. lip lead to jj,J m > so 
that /i and J m are of the same sign. The bifurcation criterion (14. ip yields 

/(Ai, X 2 )(J m -+-3-AD- /(A 2 , Ai)(A? + 2AiA 2 - \*) = 0, (4.14) 

where f(x,y) := x 3 + x 2 y + 3xy 2 — y 3 . Note that as J m — > 00 this equality 
tends to /(Ai,A2) = 0, which is the universal bifurcation criterion for neo- 
Hookean (and Mooney-Rivlin) materials [17J. Using experimental data for 
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the thoracic aorta of a 21-year old male and of a 70- year old male, Horgan and 
Saccomandi [2U] found that for these samples, J m = 2.289, 0.422, respectively. 
For these values of J m and for the three special prestrains where A3 = 1, or 
A2 = A3, or Ai = A3, we find that there exists no value of Ai such that the 
bifurcation criterion is satisfied, indicating that the aorta of the two males 
is always stable near the surface with respect to finite compressions of these 
types. Keep in mind however that the range of possible ratios is limited by 
the inequality fl4.13[) ?. For example, in plane strain A3 = 1 (no extension in 
the X3 direction), we find that 



1 + Jmil ~ V 1 + ^m 1 )/ 2 < K2 < V 1 + JmO- + V 1 + 4 -C)/ 2 - (4.15) 

For the younger aorta, this range is [0.497,2.010]; for the older, stiffer, aorta, 
the range is [0.727,1.376]. 
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